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Diversity of Chrysomelidae (Coleoptera) at a mountain range 
in the limit of the Eurosiberian region, northwest Spain: 
species richness and beta diversity 


Andrés Baselga & Francisco Novoa 


1. Introduction 


Baselga, A. & Novoa, F. 2007: Diversity of Chrysomelidae (Coleoptera) at a 
mountain range in the limit of the Eurosiberian region, northwest Spain: species 
richness and beta diversity. — Entomol. Fennica 18: 65-73. 


Chrysomelidae from the Sierra de Queixa mountains (Galicia, northwest Spain) 
were sampled, reporting 93 species. The estimated local species richness using 
several non-parametric estimators and accumulation models varies between 104 
and 142 species. To compare the Chrysomelidae fauna from Sierra de Queixa 
with other Galician assemblages we have assessed beta diversity among invento- 
ries and we have tested the differences on the zoogeographic compositions 
among areas. Sierra de Queixa is grouped with other Galician mountain ranges 
located in the transition zone between Eurosiberian and Mediterranean regions, 
and it is characterised by a high proportion of Iberian endemic species, signifi- 
cantly higher than expected. Therefore, Chrysomelidae fauna from Sierra de 
Queixa represents a unique assemblage of Eurosiberian elements that reached the 
area due to the climatic conditions and Iberian endemic species that are present in 
the Iberian mountains due to its historic role as southern refugia during glaci- 
ations. 
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tory (Baselga & Novoa 2006) as well as to assess 
beta diversity among local faunas, including in- 


Knowledge about Chrysomelidae (Coleoptera) 
from Galicia (northwest Spain) has been notably 
increased in the last years with the study of sev- 
eral local faunas (Baselga & Novoa 2000a, 
2000b, 2002a) and the publication of taxonomic 
papers (Baselga & Novoa 2000c, 2002b, 2002c, 
2003, 2004). The accumulation of such data have 
allowed to estimate the regional species richness 
and thus the completeness of the Galician inven- 


ventories from coastal environments to mountain 
ranges (Baselga & Novoa 2005). In this scope, 
we study here the leaf beetle fauna from Sierra de 
Queixa (Ourense, Spain) (Fig. 1), an area in- 
cluded in the Eurosiberian phytogeographic re- 
gion (Izco 1987, Rivas-Martinez 1987) but adja- 
cent to the limit of the Mediterranean region, 
searching for the unrecorded taxa that were pre- 
dicted to be present in Galicia by several estima- 
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Fig. 1. Location of the studied area in Galicia, northwest of the Iberian peninsula. — a. Galicia (shadowed) is lo- 
cated in the boundary between Eurosiberian and Mediterranean phytogeographic regions (discontinuous line, af- 
ter Rivas-Martinez 1987). — b. Location of Sierra de Queixa (Que) and other areas selected for the analysis of 
beta diversity. The squared pattern corresponds to the 10 km UTM grid (zone 29T, 100 km squares MH, NG, NH, 
NJ, PG, PH and PJ) and the continuous line is the boundary between phytogeographic regions (Izco 1987, 
Rivas-Martinez 1987). Each inventory is composed of several localities marked with the same symbol. See Ma- 


terials and Methods for further explanation. 


tion methods (Baselga & Novoa 2006) and sur- 
veying a new area in order to improve the descrip- 
tion of turnover patterns among Galician local 
faunas (Baselga & Novoa 2005). Our study was 
the first extensive one in Sierra de Queixa. 

The aims of this paper are (i) to describe the 
local species richness of this highly diverse 
group, (ii) to estimate the completeness of the in- 
ventory, assessing the accumulation curve gener- 
ated from our sampling, and (iii) to compare this 
fauna with other Galician ones, analysing beta di- 
versity among inventories and its differences in 
the proportions of chorotypes contributing to 
each local fauna. 


2. Material and methods 


Sierra de Queixa forms part of the Galician 
Macizo Central, a group of mountain ranges with 
summits up to 1,700 m located in the Spanish 
province of Ourense (Galicia). The range is in- 
cluded in the Eurosiberian phytogeographic re- 
gion (Izco 1987, Rivas-Martinez 1987) but Sierra 
de Queixa is adjacent to the temperate valleys of 
Sil and Bibei rivers, both belonging to the Medi- 


terranean region. Therefore, in this area there is a 
high turnover, with Eurosiberian vegetation re- 
placing Mediterranean one in a small extent. 

An extensive sampling in Sierra de Queixa 
was carried out by authors in 2002 and 2003, 
though some previous captures are also consid- 
ered in this paper. Twenty localities were sur- 
veyed (Fig. 1). They were selected in order to 
range the entire altitudinal gradient and the main 
vegetation units. A total of 789 specimens were 
collected by sweeping and beating the vegetation. 
Samplings were not standardized, and we col- 
lected all the species detected in each locality but 
not all the specimens of each species. All the ma- 
terial is deposited in A. Baselga collection in the 
Departamento de Biologia Animal, Universidad 
de Santiago de Compostela, Spain. 

To assess the completeness of the inventory 
both asymptotic models and non-parametric esti- 
mators (ICE, Chao 2 and jackknife of first order 
and second) were used (Colwell & Coddington 
1994). We used only incidence-based estimators 
since not all the specimens sawn were collected in 
each sampling and thus no abundance measures 
are available. Among them, we selected those 
recommended be Hortal et al. (2006). The accu- 
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mulation curves and the non-parametric estima- 
tors were generated with EstimateS 7.0 software 
(Colwell 2004), randomising the sample order 
100 times. Database records (Soberón et al. 2000, 
Hortal et al. 2001, Martin-Piera & Lobo 2003) 
and individuals were used as a sampling-effort 
surrogate. Our database comprised 270 records 
for 789 specimens. Each record is composed of 
the following fields: species name, locality, date, 
number of specimens and collector. Any differ- 
ence in any database field value gives rise to a 
new database record, so increments of the number 
of records provide correlative increments of the 
sampling effort (Martin-Piera & Lobo 2003). 
Thereafter, the asymptotic Clench function was 
fitted to the smoothed curve (Soberón & Llorente 
1993, Hortal et al. 2004, Jiménez-Valverde et al. 
2006): 

S,.=ae/(1t+be) (1) 
where S, is the number of species found per sam- 
pling effort unit e; a and b, the parameters of the 
function. These latter parameters were adjusted to 
the data of the curve by means of a Simplex and 
Quasi Newton method (StatSoft 2001). The as- 
ymptote (i.e. the predicted number of species) is 
calculated as a/b. 

For the comparison of the present inventory 
with other Galician faunas we have considered 
the previously published data (Baselga & Novoa 
2006), taking into account all the available inven- 
tories: coastal dunes and associated marshes 
(Dun), agricultural landscape near Santiago de 
Compostela (San), Atlantic mixed forest of low 
altitude in Fragas del Eume Natural Park (Eum), 
medium altitude mountain ranges of Dorsal 
Gallega (Dor) and Larouco (Lar), the temperate 
valley of the Sil river (Sil) and high mountain 
ranges of Ancares (Anc) and Eixo-Segundera 
(Seg). The present inventory is abbreviated as 
Que for tables and figures. Simpson’s index (B n) 
was used as beta diversity measure (Koleff et al. 
2003): 


B,,,= min(h, c) / min(b, c) +a 2) 
where a is the total number of species that occur 


in both areas, b is the total number of species that 
occur in the first area but not in the second, and c 


is the total number of species that occur in the sec- 
ond area but not in the first. This index was se- 
lected because it focuses on compositional differ- 
ences independent of species richness gradients, 
and thus it is a ‘narrow sense’ measure of turnover 
(Koleff et al. 2003). A distance matrix was con- 
structed with B „ for each pair of faunas, and 
thereafter used in a a nonmetric multidimensional 
scaling (NMDS) (Legendre & Legendre 1998), 
generated with Statistica 6.0 (StatSoft 2001). 
Three dimensions were considered because the 
screeplot showed a small change of the stress val- 
ues with dimensionality increments after this 
point (Legendre & Legendre 1998) (stress = 
0.0552). Thereafter, the matrix for these three di- 
mensions versus the nine inventories was used in 
a cluster analysis, in order to recover the major 
groups of local faunas. The tree was constructed 
with Clustan Graphics 6.03 (Clustan 2003), 
weighting the three variables after the respective 
decrease of stress explained by each dimension, 
and then grouping the faunas using Euclidean dis- 
tance and single linkage. 

Biogeographic types used in this paper were 
synthesized following the chorotypes proposed 
by Vigna Taglianti et al. (1992) and then group- 
ing them into four major categories: Iberian ele- 
ments (Ibe), for species endemic from the Iberian 
Peninsula; Mediterranean elements (Med), for 
species widespread in the Mediterranean coun- 
tries; Eurosiberian elements (Eur), for species 
widespread in Europe, or Europe and the Siberian 
range; and finally wide range elements (WR), for 
species widespread in all or a great part of the Pa- 
laearctic region, and reaching parts of both 
Eurosiberian and Mediterranean areas. These 
major divisions are established in order to make 
clear the contrast between Eurosiberian (septen- 
trional) and Mediterranean (meridional) contri- 
butions to Galician fauna which is located across 
the Eurosiberian-Mediterranean boundary. The 
other two categories are neutral regarding this as- 
pect, because almost all Iberian elements are 
present in both sides of Eurosiberian-Mediterra- 
nean limit, as well as WR elements reach both re- 
gions. 

To test the statistical significance of differ- 
ences on biogeographic patterns among invento- 
ries a chi square analysis (x°) was performed. The 
contingency table was constructed with the four 
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Med, Ibe) versus the nine inventories (rows: Dun, 
San, Eum, Dor, Lar, Sil, Anc, Seg and Que). 
Probability considered to reject the null hypothe- 
sis (differences are due to random chance) is 
95%. The degrees of freedom are 1 for the cells, 3 
for the row subtotals, 8 for the column subtotals 
and 24 for the whole contingency table, and the 
critical y” values 3.8415, 7.8147, 15.5073, 
36.4150 respectively 


A total of 93 species of Chrysomelidae have been 
collected in Sierra de Queixa. The full list of spe- 
cies and exact sampling localities are available 
from the authors on request. The leaf beetle fauna 
in Sierra de Queixa is mostly composed of WR el- 
ements (43.0%). The proportions of Eur and Med 
elements are equal (16.1%), and the Ibe ende- 
misms reach near a quarter of the inventory 
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Table 1. Estimated species richness using database 
records and individuals as sampling effort surrogate. 








Estimate Records curve Individuals curve 
Clench 142 105 
ICE 115 105 
Chao 2 115 104 
Jacknife 1 124 111 
Jacknife 2 135 116 





(24.7%). Among the cited taxa, 6 species are spe- 
cially interesting: Psylliodes cervinoi Baselga & 
Novoa is firstly recorded here since its descrip- 
tion from Sierra Segundera (Baselga & Novoa 
2003), extending its known distribution to 
Queixa mountain range, whereas Stylosomus mi- 
nutissimus (Germar, 1823), Pachybrachis suffri- 
ani Schauffus, 1862, Galeruca interrupta (Illiger, 
1802), Luperus longicornis (Fabricius, 1781) and 
Psylliodes obscuroaeneus Rosenhauer, 1856 are 
new records for Galicia, increasing the regional 
inventory from 276 to 281 species. 

The completeness of the present inventory 
was evaluated by means of accumulation models 
and several non-parametric estimators, using 
both database records and individuals as sam- 
pling effort surrogates (Fig. 2). Our results show 
that a total richness between 104 and 142 species 
could be expected in the studied area, depending 
on the curve assessed and the estimator selected 
(Table 1). Therefore between 65% and 90% of the 
estimated number of species of Chrysomelidae 
living in Sierra de Queixa seems to be detected in 
our field study, but we should expect that between 
12 and 49 taxa will be added to this inventory in 
the future. Moreover, in both curves the slope at 
the end of the survey is quite low: 0.1195 using 
records and 0.0156 using individuals. These val- 
ues imply that more than 8 records or 64 individu- 
als should be added to the inventory, and thus the 
completeness of the inventory is quite high 
(Jiménez- Valverde & Hortal 2003, Hortal et al. 
2004). 

To compare the leaf beetle fauna from Sierra 
de Queixa with those of other Galician areas we 
have assessed beta diversity among inventories, 
considering the presence data reported in Baselga 
and Novoa (2006). The NMDS (Fig. 3a) was per- 
formed with three dimensions (stress=0.0552) 


g uojsuawiq 
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Fig. 3. Comparison of leaf beetles from Sierra de 
Queixa with other Galician local faunas. — a. Non- 
metric multidimensional scaling (NMDS) based on 
Simpson’s measure of beta diversity (stress = 
0.0552). — b. Dendrogram from the cluster analysis, 
constructed after the matrix for faunas vs. dimensions 
produced by the NMDS. Euclidean distance and sin- 
gle linkage were used. 


because the screeplot showed a small change of 
the stress values with dimensionality increments 
after this point (Legendre & Legendre 1998). Ma- 
jor groups of Galician faunas were recovered by a 
cluster analysis, using the matrix for faunas vs. 
dimensions produced by the NMDS. The dendro- 
gram constructed from this matrix shows three 
major faunistic groups (Fig. 3b). The leaf beetle 
fauna from Sierra de Queixa is grouped with 
those of other mountain ranges which are located 
in the trasition zone between Eurosiberian and 
Mediterranean regions (Ancares, Larouco and 
Eixo-Segundera). The other two clusters were 
also found elsewhere (Baselga & Novoa 2005): 
on one hand, the low and medium altitude areas 
closer to the Atlantic coast (coastal dunes, Santi- 
ago, Eume Natural Park and Dorsal Gallega) and, 
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Table 2. Contingency table for the biogeographic patterns versus the local inventories. Significant differences are 
noted in bold face. Observed and expected values are noted as (O) and (E), respectively. X total value for the 


whole table = 49.19. 








Corotype WR WR% Eur Eury° Med Medy? Ibe Ibex? Sub- Subt. x’ 
total 

Dun (O) 72 11 8 113 

Dun (E) 595 26 18.2 29 18.5 0.7 16.8 4.6 10.8 
San (O) 66 16 12 104 

San (E) 547 2.3 16.8 0.0 17.0 29 155 0.8 6.0 
Eum (O) 43 18 17 87 

Eum (E) 45.8 0.2 14.0 1.1 14.2 1.9 12.9 1:3 4.5 
Dor (O) 45 14 11 83 

Dor (E) 43.7 0.0 13.4 0.0 13.6 0.0 124 0.2 0.2 
Lar (O) 40 11 10 73 

Lar (E) 38.4 0.1 11.8 0.1 1.9 0.0 10.9 0.1 0.2 
Sil (O) 50 10 13 103 

Sil (E) 542 0.3 16.6 2.6 16.9 10.3 15.3 04 13.6 
Anc (O) 69 31 21 141 

Anc (E) 74.2 04 228 3.0 23.1 0.4 21.0 0.0 3.8 
Seg (O) 35 15 15 77 

Seg (E) 40.5 0.8 12.4 0.5 12.6 0.0 11.5 1.1 2.4 
Que (O) 40 15 23 93 

Que (E) 49.0 1.6 15.0 0.0 15.2 0.0 13.8 6.1 77 
Subt. (O) 460 8.31 141 10.26 143 16.2 130 14.4 874 49.2 





on the other, the fauna from the Sil Valley that is 
the most different one. 

The comparison of the biogeographic compo- 
sitions of the local inventories shows similar re- 
sults (Table 2). The column subtotal y° for Sierra 
de Queixa is quite high (close to the critical value 
but not significant at p<0.05) almost completely 
due to the large proportion of Iberian endemics, 
which is significantly higher than expected by 
chance. Only two inventories are significantly 
different to the others: Sil valley, due especially to 
its high proportion of Mediterranean elements 
and coastal dunes due mainly to its small percent- 
age of Iberian endemisms. These disparities pro- 
duce a high total y? value for the whole table, im- 
plying significant differences among the com- 
pared faunas. 


4. Discussion 


Chrysomelidae constitute an interesting group to 
develop biodiversity estimates due to two main 
reasons: on the one hand it is a hyperdiverse fam- 
ily and, on other, they provide a good representa- 


tion of the free living phytophagous functional 
group within arthropods. The study of leaf beetles 
from Sierra de Queixa has revealed a very inter- 
esting local fauna, quite singular among the re- 
gional assemblages due to the striking high pro- 
portion of Iberian endemics (24.7%) (Fig. 4, Ta- 
ble 2) or the citation of several taxa new for 
Galicia. The addition of five new species to the 
regional inventory may be considered valuable, 
since completeness of Galician Chrysomelidae 
list was already very high, between 85 and 95% 
of the estimated richness (Baselga & Novoa 
2006). This saturation is more difficult to get at a 
smaller scale because of the great number of 
samplings that would be required for each local 
area, and thus a degree of completeness between 
65 and 90% may be considered acceptable, con- 
sidering that we are working with a highly diverse 
group of invertebrates. Therefore, further sam- 
pling effort in Sierra de Queixa would add new 
taxa to the inventory (between 12 and 49 taxa af- 
ter our estimates) but current knowledge of 
Queixa leaf beetles is already suitable in order to 
describe the local assemblage and the regional 
biodiversity patterns. Comparing the different es- 
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timates (Table 1), the maximum value is recov- 
ered adjusting the Clench function to the accumu- 
lation curve using database records as sampling 
effort surrogate (142), whereas the minimum one 
is provided by the Chao 2 estimator in the individ- 
uals curve (104). In all cases, estimates of the 
curve based on individuals are lower than those 
yielded by the curve based on records, but the dif- 
ference between the Clench asymptotes of both 
curves is the most remarkable discrepancy. Ex- 
amination of the frequency of records per species 
in the curve based on records shows an important 
number of species with 4 or less records (not 
shown). This causes the curve to overestimate to- 
tal species number. Such drawback is less impor- 
tant when using individuals. In this curve, values 
of Clench asymptote and non-parametric estima- 
tors are more similar. More research on the be- 
haviour of species accumulation curves is 
needed, because in our data set the number of sin- 
gletons and doubletons is considerable, so non- 
parametric estimators should be overestimating 
the richness value, and thus it is striking that as- 
ymptotic model yield grater values, contrary to 
the results found by Hortal et al. (2006). 

When the Chrysomelidae fauna from Sierra 
de Queixa is compared with other Galician inven- 
tories analysing beta diversity and ordering local- 
ities by means ofa NMDS, three major groups are 
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found. Queixa is grouped together with other 
mountain ranges located in both sides of the 
Eurosiberian-Mediterranean phytogeographic 
boundary (Izco 1987, Rivas-Martinez 1987). 
This group was previously recovered in similar 
comparison performed without including present 
data (Baselga & Novoa 2005) and the addition of 
Queixa inventory to this cluster is not surprising 
taking into account its location. The conclusion 
that faunistic isolation of Sil Valley (Figs. 3 and 4) 
is due to environmental factors (i. e. temperature, 
pluviosity), because no geographical barriers are 
interposed, is here newly supported. Sierra de 
Queixa and Sil Valley are strictly contiguous as 
well as Sil Valley and Eixo-Segundera. Also the 
clustering of two faunas located in the phyto- 
geographic Mediterranean region (Larouco and 
Eixo-Segundera) with the Eurosiberian invento- 
ries is again corroborated. Hence, it should be 
emphasized that although the Larouco and Eixo- 
Segundera ranges are included into the phyto- 
geographic Mediterranean region together with 
the Sil Valley, only the latter area has an unambig- 
uous Mediterranean-type Chrysomelidae fauna 
(Fig. 5). This conclusion was pointed in Baselga 
& Novoa (2005) but it is statistically tested here 
for the first time. The contingency table shows 
significant differences in the distributions of bio- 
geographic categories among local faunas (Table 


72 Baselga & Novoa * ENTOMOL. FENNICA Vol. 18 


2). In the case of Sil valley Mediterranean ele- 
ments are significantly higher than expected. 

The zoogeographical zonation is different for 
each group considered, but as a general pattern 
Eurosiberian faunas occupy areas that Rivas- 
Martinez (1987) includes in the Mediterranean 
region in his classification of Iberian vegetation, a 
pattern also pointed out by Hermida et al. (1994) 
for western Iberian gastropods, Rodriguez et al. 
(1997) for western Iberian earthworms and 
Hortal-Mufioz et al. (2000) for Portuguese dung 
beetles. As Williams (1996) comments, the arbi- 
trary “border lines” between biogeographic re- 
gions are not realistic representations of the eco- 
logical reality, because they can be wider than the 
grain used for their delimitation. In fact, it seems 
that more than a sharp boundary, the transitional 
zone is composed by a mosaic of leaf beetle fau- 
nas depending on the local climatic conditions. 
These environmental variables are mostly deter- 
mined by the complicated regional geomorphol- 
ogy. Therefore, mountains located south of the 
boundary (and then included in the Mediterra- 
nean region; 1.e. Larouco, Eixo-Segundera) host 
Eurosiberian faunas, whereas temperate valleys 
located north to the boundary shelter Mediterra- 
nean faunas (i.e the Navia valley in Ancares 
range; Baselga & Novoa 2000b). This mosaic 
produces strong gradients and sharp faunal re- 
placements in some places. This is the case of Si- 
erra de Queixa and the contiguous Sil Valley, 
where the transition from unambiguous Euro- 
siberian localities in the mountains (close to 1700 
masl) to clearly Mediterranean areas in the tem- 
perate valleys (at 300 masl) occurs in less than 10 
km. Considering that Mediterranean and Iberian 
endemic species could be relict taxa isolated in 
Pleistocene refugia during glaciations (see Ribera 
& Vogler 2004), the faunistic uniqueness of Si- 
erra de Queixa is based in the confluence of both 
Eurosiberian taxa (able to attain the mountain 
range due to current climatic conditions) and rel- 
ict Iberian species (which are present in the 
mountain range due to its historic role as southern 
refuge during glaciations). 

In summary, the check-list of Chrysomelidae 
from Sierra de Queixa has reported interesting 
faunistic findings as the new record of Psylliodes 
cervinoi, the citations of new taxa for Galicia or 
the description of this Iberian endemic-rich 


fauna, unique among regional assemblages. 
Moreover, the addition of this study to the re- 
gional knowledge has allowed us to reassess leaf 
beetle beta diversity in the region. The new data 
support the conclusions reported in our study of 
the faunistic heterogeneity of Galician leaf bee- 
tles (Baselga & Novoa 2005) and it should be as- 
sumed that the described patterns are robust. 
They show that zoogeographic distribution of 
Chrysomelidae faunas in northwestern Spain do 
not fit exactly the standard phytogeographic 
zonation (Izco 1987, Rivas-Martinez 1987). Only 
Sil Valley has an unambiguous Mediterranean- 
type fauna, whereas mountain ranges of Larouco 
and Eixo-Segundera, though located in the 
phytogeographic Mediterranean region, are 
grouped with Eurosiberian leaf beetle faunas. 
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